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A recent paper by Zeidler et al. |l| describes a neu- 
tron scattering experiment on water in which oxygen 
isotope substitution is successfully achieved for the first 
time. Differences between scattering patterns with dif- 
ferent oxygen isotopes give a combination of the 0-0 
and O-H (or O-D) structure factors, and the method ele- 
gantly minimises some of the problematic inelasticity ef- 
fects associated with neutron scattering from hydrogen. 
Particular conclusions of the new work are that the OH 
bond length in the H2O molecule is about 0.005A longer 
than the same bond in D 2 0, and that the hydrogen bond 
peaks in both liquids are at about the same position. 

Notwithstanding the substantial progress demon- 
strated by the new work, the comparison with our own 
results |2J by Zeidler et al. is in our opinion misleading. 
The last paragraph of [l| states that the new experimen- 
tal results "...originate directly from measured data sets 
and not from a refinement of models using the diffrac- 
tion data, wherein different starting points can lead to 
quite different conclusions [36,41]" Firstly the results 
reported in references 36 and 41 Q do indeed originate 
directly from measured datasets, as is abundantly clear 
from the papers themselves. Equally the reported OH 
bond lengths in Q are a direct consequence of the par- 
ticular diffraction data presented, as highlighted in Fig. 2 
of that paper, and do not depend on the computer simu- 
lation method used. Both papers use different empirical 
potential structure refinement (EPSR) computer simu- 
lations to fit radiation scattering data and derive 00, 
OH and HH radial distribution functions compatible with 
those data. Comparing Fig 3. of [1] with Fig. 6 of Q and 
Fig. 6 of [5], which include independent datasets from 
both x-ray and neutron (reactor and pulsed) sources, and 
from "different starting points", the radial distribution 
functions for water extracted by this method from all 
these different datasets are in fact very similar. But they 
are of course not identical and this is known to be due to 
systematic uncertainties which are difficult to quantify, 
particularly where scattering from hydrogen or small iso- 
tope differences is involved [5|. Zeidler et al. distance 
themselves from the EPSR approach, yet by choosing 
one potential, TTM3-F, out of several possibilities they 
have in fact performed what EPSR does automatically, 
namely selected a potential that reproduces their data. 
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FIG. 1: Simulated first order difference functions (solid 
lines), AF D (Q) = 0.0059[S OD (Q) - 1] + 0.00262[5 OO (Q) - 
l](top), AF H (Q) = -0.0033[S O h(Q)-1] + 0.00263[S , oo(Q) 
1] (bottom), based on the EPSR simulations reported in 
The oxygen isotope substitution data from [l| are shown as 
open circles. The inset shows the high Q region in more detail. 
Curves are shifted vertically for clarity. 



This simulation is then used to claim that the hydrogen 
bond length in H2O is about the same length as that in 
D2O, just as in [2j the EPSR simulation is used to claim 
these bonds are different. Note however that the results 
reproduce the change in x-ray scattering pattern be- 
tween heavy and light water [f|, something that TTM 
potentials tend to underestimate 0j- 

To illustrate the ambiguities associated with using any 
particular dataset to characterise the structure of wa- 
ter we have constructed the oxygen first order difference 
functions, AF H (Q) and AF D (Q) (Figure 1(a) of Q) us- 
ing the simulated 00 and OH partial structure factors 
that were extracted in the previous analysis [2j , weighted 
according to the new values of the oxygen isotope scat- 
tering lengths, and without refinement against the new 
data (Figure [1} . Although there is some disagreement 
in the region of the main peak near Q = 2A _1 in the 
AFd(Q) function, it is striking how close the earlier 
EPSR simulations approach the new data. A similar dis- 
crepancy in AFd(Q) between simulation and experiment 
near Q = 2A _1 is observed in Figure 1(a) of [l|. 
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